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When a resonance associated with electromagnetically in- 
duced transparency (EIT) in an atomic ensemble is modulated 
by an off-resonant standing light wave, a band of frequencies 
can appear for which light propagation is forbidden. We show 
that dynamic control of such a bandgap can be used to co- 
herently convert a propagating light pulse into a stationary 
excitation with non-vanishing photonic component. This can 
be accomplished with high efficiency and negligble noise even 
at a level of few-photon quantum fields thereby facilitating 
possible applications in quantum nonlinear optics and quan- 
tum information. 

PACS numbers 03.67.-a, 42.50.-p, 42.50.Gy 



Techniques for coherent control of light-matter inter- 
action are now actively explored for storing and manupi- 
lating quantum states of photons. In particular, using 
electromagnetically induced transparency (EIT) p]J| and 
adiabatic following of "dark-state polaritons" [BL the 
group velocity of light pulses can be dramatically de- 
celerated and their quantum state can be mapped onto 
metastable collective states of atomic ensembles 0. 

In contrast to such a coherent absorption process, the 
present Letter describes how a propagating light pulse 
can be converted into a stationary excitation with non- 
vanishing photonic component. This is accomplished via 
controlled modification of the photonic density of states 
in EIT media by modulating the refractive index with an 
off-resonant standing light wave. By varying the prop- 
erties of the resulting photonic band structure in time, 
the original light pulse can be converted into an exci- 
tation inside the bandgap where its propagation is for- 
bidden. Long storage of excitations with non-vanishing 
photonic component may open interesting prospects for 
enhancement of nonlinear optical interactions |^||. In 
particlular, an intriguing and practically important 
application of this effect for interactions between few- 
photon fields is dicussed in the concluding paragraph of 
this Letter. 

Before proceeding, we note that there exists a sub- 
stantial litterature on photonic bandgap [pj materials. 
Recently photonic bandgap structures have been inves- 
tigated theoretically |10| for strong coupling of single 
atoms with photons. Photonic bandgap based on inter- 
action with atoms in an optical lattice were also investi- 
gated . We also note other related work on EIT-based 
control of the propagation properties of light in atomic 
media @. 

The key idea of the present approach can be qualita- 



tively understood by first considering a medium consist- 
ing of stationary atoms with a level structure shown in 
Fig. la. The atoms are interacting with a weak signal 
field and two strong fields. The running wave control field 
£! c is tuned to resonant frequency of the \b) — > |c) tran- 
sition. In the absence of the field f2 s , this situation cor- 
responds to the usual EIT: in the vicinity of a frequency 
corresponding to two-photon resonance the medium be- 
comes transparent for a signal field. This transparency is 
accompanied by a steep variation of the refractive index. 

The dispresion relation can be further manipulated by 
applying an off-resonant standing wave field with Rabi 
frequency £l s {z) — 2£l s cos(k s z) and a frequency detun- 
ing A. This field induces an effective shift of resonant 
frequency (light shift) that varies periodically in space, 
resulting in a spatial modulation of the index of refrac- 
tion according to 5n(z) = (c/v g )A-^ cos 2 (k s z), where 
c/Vg is the ratio of speed of light in vacuum to group 
velocity in the medium. When the modulation depth 
is sufficiently large, signal light propagating near atomic 
resonance in the forward z direction with wavenumbcr 
k near k s may undergo Bragg scattering into the back- 
ward propagating mode with wavenumber —k. In direct 
analogy to e.g., optical interferometers, the scattering of 
the counterpropagating fields into each other can modify 
the photonic density of states. In particular, a range of 
frequencies ("photonic bandgap") can appear for which 
light propagation is forbidden jl3| . According to a stan- 
dard technique to analyze the resulting band structure, 
Bloch's theorem can be applied so that the propagating 



solutions obey E(z+a) 



a E(z), where K is the Bloch 



wave vector. Imposing this condition and assuming that 
the wave vectors of the fields are close (k ~ fc s ), we can 
solve for the band structure and obtain near two-photon 
resonance 



cos(Ka) = cosh 



^«VA^-(,-*) 2 ) , (1) 



where g = p 



2 h " y is the atom-field coupling constant, 

N is the number of atoms, A s = fi^/A is the amplitude 
of the light shift modulation, p is the dipole moment of 
the a — b transition, V the quantization volume and the 
factor g 2 N/fl 2 corresponds to c/v g . For frequencies such 
that \u> — u>b a \ < |A S | a bandgap is created: the Bloch 
wavevector acquires an imaginary part and the propaga- 
tion of waves in the medium is forbidden. For an outside 
observer such a medium can be viewed as a mirror: an 
incident wave with frequency inside the bandgap would 
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undergo almost perfect reflection. Calculations in Fig. 2 
indicate that this qualitative result remains valid even for 
realistic EIT conditions, including a finite transparency 
bandwidth and finite ground-state decoherence rate jbc- 

A specific, distinguishing feature of the present scheme 
is the possibility of dynamically changing the properties 
of the medium by switching in time the fields £l s (t) and 
Q c (i) on and off. In particular, by combining the tech- 
niques of Q with the present idea, the following scenario 
can be implemented: first, with the standing wave turned 
off Q s = 0, a forward propagating pulse is stored in the 
medium as a Raman coherence between levels \b) and 
|c). Then, by switching on both the "control" field tt c 
and the standing wave field fl s , the pulse can be released 
into the bandgap and initially propagates in the forward 
direction. In the presence of a bandgap, the forward 
(£+) and backward {£-) components are coupled due to 
Bragg scattering off the index grating, so that amplitude 
in the forward mode +k is converted into amplitude in 
the backward mode —k and vice versa. In this case, the 
pulse can be effectively trapped in the photonic bandgap. 

We now turn to a detailed description of the dynamic 
trapping procedure. We are interested here in the prop- 
agation of fields with possibly non-trivial statistics, such 
as single photon fields, so that a quantum description is 
used. In the presence of the standing wave it is conve- 
nient to decompose the propagating signal fields into two 
slowly varying components £+ (propagating forward) and 
£_ (propagating backward) so that the electric field is 

= ^v\^a=±^^y {v/c)(a ^ ct) +h-c] and 
couples resonantly to the transition between the ground 
state | a) and the excited state \b), the carrier frequency 
of the optical field being v = u) a f,. Two time dependent 
classical driving fields with Rabi-frequencies Q s (t) and 
fl c (t) are used to control the propagation as shown if 
Fig. la. 

To describe the quantum properties of the medium, 
we use collective slowly varying atomic operators fli] ] 

o>„(z,t) = Y,j=i \N)( v i\ er ' lu ''"' t where the sum is 
performed over a small but macroscopic volume contain- 
ing Nz 3> 1 atoms around position z. The interaction 
hamiltonian is then, in a rotating frame 



H 



Aa dd 



+ n c e lkcZ a ac + 2Q S cos(k s z)a dc + h.c] } , 



(2) 



where fco —vjc and N the number of atoms. 

Since the two propagating fields £+ and £- interact 
with the atoms, we expect an optical coherence a\, a to 
appear as well as a Raman coherence <7f, c . Moreover 
these two fields will give rise to coherences with distinct 
spatial variations, i.e., varying as e ik ° z for the compo- 
nent of &ba induced by £+ while that due to £- will vary 
as e~ lk ° z . We therefore decompose the optical and Ra- 
man coherences according to these two distinct spatial 



variations &ba(z,t) = a^ a (z,t)e tk ° z +a ba (z,t)e tk ° z and 
a bc (z,t) = a+J^z,ty( k °- k -> +a~ c (z,t)e- l( - ko+k ^ z . Using 
slowly varying envelopes, we then have the equations of 
motion for the forward and backward modes 



Q- t ±c—)£±(z,t) = igNa b i a (z,t). 



(3) 



Assuming weak quantum fields and solving perturba- 
tively, we find to lowest order in the weak fields and in 
an adiabatic approximation (assuming Q s (t) and O c (i) 
change in time slowly enough ||) 
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(4) 
(5) 



where Ak = k s — k , A s = |£! S | 2 /A is the amplitude of 
the spatially modulated light shift caused by the stand- 
ing wave field n s and F±(t) are (^-correlated noise forces. 
Note that in the adiabatic limit the noise forces are neg- 
ligible H . The propagation equations are thus 



1*4 1 f '■•■" = 



S 2 N B £ ± {z,t) 
fL dt fl c 



+ j 9%N a / t ^ Z '^ c ±2iAkz 



a. 



fir 



(6) 



which indicates that the forward and backward slowly 
propagating modes become coupled. Specifically, the first 
term on the right-hand side gives rise to propagation at 
the group velocity v g = c/(l + ^?) fl5| while the sec- 
ond term gives rise to coupling between the forward and 
backward propagating modes. This coupling is optimum 
when the effective phase matching Afc = k s — fco = is 
achieved jl3| . Note that both the "control" field 17 c and 
the standing wave amplitude f2 s can be time-dependent 
and that as long as changes are slow enough (adiabatic 
limit H) the above equations describe the correct dy- 
namics of the coupled modes. 

To obtain a solution in the case of time-dependent 
fields r2 c (t) and 0, s (t), we introduce new quantum fields 
i& + (z t t) and ^>^(z,t) (forward and backward propagat- 
ing dark-state polaritons ||) ^±(z,t) — cos 6{t)£± (z, t) — 
sin 6 '(t)-\/Nafa,(z,t), where tan 2 6{t) = is the mix- 

ing angle between the photon and matter components 
of the polariton. The polaritons then obey the coupled 
equations 



d_ ± d_ 

dr dz 



*± = iA s tan 2 6>(i)* T , 



(7) 



where r(t) = f dt' cos 2 9(t'). Eq. (^) describes prop- 
agation with velocity v g (t) — c cos 2 6(t) of the two po- 
laritons (traveling in opposite directions) and coupling 
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with rate A s (t) sin 2 9{t). Note that in the limit c 3> v g , 
the photonic component £± ~ (Q c /g\/N)'&± is finite for 
non-zero control field f2 c . 

We consider now the scenario in which the standing 
wave beams are initially off and the control field is on, 
with Rabi frequency f2™ (corresponding to a group ve- 
locity v g n ). A forward propagating photon wavepacket 
can then be stored in the medium in the form of a Ra- 
man coherence <7^,(z,i) and subsequently released jjj. 
We consider the case when the standing wave field is 
first switched on, establishing the bandgap, followed by 
the control field (with Rabi frequency correspond- 
ing to a group velocity v g , possibly different from v l g n ), 
releasing the pulse in the bandgap medium. For sim- 
plicity we consider the case when the standing wave is 
switched on before or simultaneously with the control 
field, so that, the coupling rate A s (r) tan 2 9{t) does not 
depend on r. In this case, we solve ([?]) by Fourier trans- 
forming ^±(z) = 7^- j dk e lkz i&±(k) to obtain 



^ + (fc, t) = cos(£t) — i-^- sin(^r) 
$_(fc,r)=-i^Hin(CT)fc + (M), 



(8) 



where \ = AJt) tan 2 9{t) and £ = \J (he)* + 
According to (|8|), the various Fourier components 
(wavenumber k) of the pulse cycle back and forth be- 
tween the corresponding forward and backward modes 
at a rate which depends weakly on k. In particular when 
the spatial extent of the pulse inside the medium is large 
enough to give a relevant range of wavenumber negligible 
compared to the strength of the coupling of forward and 
backward modes, pulse distortion is negligible and the 
spatial envelopes have the time dependence ^ + (z,t) — 
cos(xr)# + (z,0) and#_(z,r) = -i sin(xr)#+(z, 0). The 
wavepacket periodically cycles between a forward and 
backward propagating component, the result of which 
is trapping of the pulse in the medium as shown in Fig. 
3. The wavepacket is trapped as a combination of light 
pulse and Raman coherence. 

The above analysis involves an adiabatic approxima- 
tion and ignores the decay of Raman coherence. In order 
to ignore the motion compared to the coupling we require 
that j;>fc and since the maximum k can be estimated 
from the initial length of the pulse in the medium we 

find that this is equivalent to requiring that A S T 3> 
where T is the duration of the initial pulse. As seen 
from (|J) the effect of non-zero values of k is that the 
trapped pulse become spatially distorted. Expanding 
Vx 2 + (kc) 2 r ~ [ X + (fec) 2 /(2x)]r we need r < x/{kc) 2 7 
which gives after expressing r in terms of real time t that 
Hfr- <C A s T(v g n /v g ) 2 , where t int is the maximum time 
during which the pulse may be trapped without suffering 
distortion. Furthermore, taking into account the limits 
imposed by adiabaticity (i.e., modulation of index occurs 



within the transparency window A s <§; (fl®) 2 /j) and the 
fact that the trapped pulse must fit inside the medium 
when travelling at the reduced group velocity, we find 
that the trapping or interaction time is limited by 



< 



9 2 N v? 1 ^ 
7 a6 (c/L) v° g ' 7 fcc ' 



(9) 



where L is the length of the medium. This limiting quan- 
tity corresponds to the density length-product and can 
be rather large for optically dense medium. 

To summarize, we have shown that by spatially mod- 
ulating the dispersive feature of the EIT resonance it is 
possible to induce a photonic bandgap. By dynamically 
controlling the resulting band structure, a propagating 
light pulse can be converted into a stationary excitation 
which is effectively trapped in the medium. 

To conclude, we note some interesting avenues opened 
by this work. First, we note that the present work is 
not restricted to the use of stationary or cold atoms, for 
example a Doppler-free configuration involving pairs of 
copropagating fields is shown in Fig. lb. In this case, the 
two polaritons are asociated with distinct atomic states 
|c) and |c'). Each polariton corresponds to a Doppler-free 
Raman configuration and they are coupled by a Doppler- 
free two-photon transition. Second, this work may open 
interesting prospects for nonlinear optics. For example, 
a trapped photonic excitation can be used to induce a 
light shift via interaction with another atom-like polari- 
ton. Large nonlinear phase shifts at the single-photon 
level can be expected and open up the way for possi- 
ble applications in quantum non-linear optics and quan- 
tum information without the limitations associated with 
traveling wave configurations |6[|l6}| and without invoking 
cavity QED techniques []l7f . Finally, it is intriguing to 
consider extension of these ideas to manipulate photonic 
bandgap in condensed matter. 

This work was supported by the NSF through ITR 
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FIG. 1. Atomic level configuration for EIT-induced pho- 
tonic bandgap. a) Stationary atoms scheme, b) moving atoms 
scheme. The standing wave of Rabi frequency Q a is detuned 
by A from resonance with the |c) — * \d) transition, giving rise 
to a spatially modulated light shift A B (z,t) — \Q s (z,t)\ 2 / A. 
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FIG. 2. a) Dispersion relation oj(K) for waves propagating 
in EIT-induced bandgap medium. Plotted is the detuning 
from two-photon resonance uj = v — vq in units of linewidth 7 
vs. the Bloch wavevector K in units of 7/c. For \u\ 3> 7, the 
dispersion relation corresponds to the usual v = cK while 
near resonance the slope is changed to the group velocity 
v — v g K. For < A s a bandgap appears, as shown in the 
inset (full line: real part of K — k s , dashed line: imaginary 
part), b) Reflection coefficient vs. frequency near EIT reso- 
nance for a sample of 87 Rb atoms corresponding to a density 
n ~ 10 13 cm -3 and length of medium L ~ 4cm (jbc = 1kHz, 
7 ~ 10MHz, gy/N ~ 400MHz, c/v g ~ 10 3 and A s ~ 400kHz), 
giving a peak reflectivity of 98.9%. 
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FIG. 3. a) Amplitude of forward and backward propagat- 
ing polaritons 9+(z,t) and ty-(z,t). b) Corresponding elec- 
tric fields and c) total intensity (forward and backward com- 
ponents) averaged over the optical wavelength. Also shown 
is the time-dependence of the "control" field fl c (t) (dotted 
line) and of standing wave field Q s (t) (dashed line). Note 
that v l g n /Vg ~ 15 here so that initial motion of the pulse is 
noticeable on these plots. Axes are in arbitrary units. 
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